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Abstract : 

axisymmetric  body,  designed  for  transition  research  in  a  1.22m  diameter 

water  tunneL  is  described.  Boundary  layer  transition  data  are  presented 

as  functions  of  the  heating  power  supplied  to  the  body  and  the  total 

concentration  of  free-stream  particulate  matter  in  the  water.  Body  surface 
temperatures  range  from  0  to  25°C  over  the  ambient  water  temperature,  and 
the  total  heat  supplied  ranges  from  0  to  93.3  kW.  Transition  arc  length 
Revnolds  numbers  are  found  to  varv  from  A. 3  million  for  the  bodv  operating 
coLd  to  36.4  million  for  the  maximum  heat  level  considered.  The 
concentration  of  free-stream  particles  is  shown  to  affect  transition 
RevnoLds  number.  These  particles  range  in  diameter  from  10  to  70  pm  and 
tneir  concentration  ranges  from  less  than  5  to  198  particles  per  cm3. 

The  decrease  in  transition  Revnolds  number  due  to  the  higher  concentration 
of  particles  is  of  order  30  percent. 
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Early  wind  tunnel  experiments  such  as  those  of  Frick  and  McCullough  (1942) 
and  Liepmann  and  Fila  (1947)  suggest  that  the  transition  location  of  a  flat 
plate  boundary  layer  in  air  is  moved  forward  as  a  result  of  surface  heating. 
Because  this  effect  is  suspected  to  be  due  to  the  variation  of  fluid  viscosity 
with  temperature,  it  has  long  been  speculated  that  the  opposite  effect  would 
occur  in  water.  In  particular,  the  viscosity  of  water  decreases  with 
increasing  temperature  which  results  in  a  viscosity  gradient  close  to  a 
heated  wall.  The  lower  viscosity  near  the  wall  would  tend  to  promote  larger 
velocities  near  the  wall  which  results  in  a  fuller,  more  stable  laminar 
velocity  profile  than  the  Blasius  profile.  Cooling  in  water  (and  heating  in 
air)  results  in  higher  viscosity  near  the  wall  which  can  lead  to  an  inflected 
velocity  profile,  which  is  unstable. 

Numerical  calculations  by  Wazzan,  Okamura,  and  Smith  (1968,  1970)  confirm 
that  wall  heating  can  produce  large  increases  in  the  transition  Reynolds 
number  of  water  boundary  layers.  In  this  analysis,  the  stability 
characteristics  of  the  boundary  layer  is  predicted  from  the  solution  of  a 
modified,  fourth-order  Orr-Sommerf eld  equation.  Lowell  and  Reshotko  (1974) 
refined  these  predictions  further  by  introducing  a  coupled  sixth-order  system 
of  vorticity  and  energy  disturbance  equations.  Such  calculations  have 
demonstrated  that  transition  Reynolds  numbers  (based  of  free-stream  velocity 
and  distance  from  the  leading  edge  of  a  flat  plate)  can  reach  as  high  as 
200  million  for  wall  temperatures  that  are  40°C  above  the  water  temperature. 
These  analytical  methods  and  others  are  thoroughly  reviewed  by  Reshotko 


(  1976) . 
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From  a  practical  point  of  view,  the  delaying  of  transition  can  result  in 
significant  drag  reduction.  Reshotko  (1978)  presents  a  series  of  calculations 
of  skin  friction  reductions  for  a  heated  flat  plate  at  zero  incidence  .  Tnese 
calculations  account  for  the  effectiveness  of  transmitting  the  reject  heat 
from  a  propulsion  system  to  the  water  in  the  desired  manner.  Even  with  a 
moderately  efficient  heat  exchanger,  drag  reductions  in  excess  of  60  percent 
are  possible  at  transition  Reynolds  numbers  of  the  order  of  70  million. 

The  potential  payoff  in  drag  reduction  has  triggered  a  fair  amount  of 
recent  experimental  research  in  heated  laminar  boundary  layer  flows. 

Strazisar,  Reshotko,  and  Prahl  (1977)  present  experimental  results  for  the 
stability  of  an  artificially  excited  laminar  layer  over  a  heated  flat  plate  at 
zero  incidence.  For  low  wall  overheats  (AT  =  Tw  -  T0  =  2.8°C,  where 
Tw  =  wall  temperature  and  T0  =  water  temperature),  they  establish  neutral 
stability  curves  from  which  critical  Reynolds  numbers  can  be  determined. 

The  results  are  in  reasonable  agreement  with  predictions.  This  work 
represents  one  of  the  first  verifications  of  the  predictions  of  Wazzan,  et  al 
t  968)  and  Lowell  and  Reshotko  (1974). 

The  experiments  of  Strazisar,  et  al  (1977)  are  limited  to  very  low 
streamwise  Reynolds  numbers  and  overheats.  Barker  and  Gile  (1981)  describe 
a  series  of  heated  pipe  experiments  of  which  one  of  the  goals  is  to  explore 
the  upper  limits  of  transition  Reynolds  number  through  wall  heating.  Their 
experiments  are  conducted  on  the  inside  of  a  long  (6.1m),  electrically  heated 
pipe.  Because  the  displacement  thickness  is  thin  relative  to  the  pipe  radius, 
the  boundarv  Laver  development  is  approximately  the  same  as  that  of  a  zero- 
pressure  gradient  flat  plate.  For  a  walL  overheat  of  8°C,  Barker  and  Gile 
measure  a  transition  Revnolds  number  of  47  miLLion.  As  the  wall  heating  is 
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Increased  further,  these  authors  note  that  there  is  no  further  improvement  in 
transition  Reynolds  number,  which  is  in  contradiction  with  theory.  Barker  and 
Gile  attribute  the  observed  limitations  to  flow  asymmetries  and  free-stream 
particulates.  The  flow  asymmetries  originate  in  the  entrance  region  of  the 
pipe  and  are  in  the  form  of  a  large  longitudinal  vortex  pair.  Yao  (1977) 
analyzes  this  flow  in  the  context  of  the  flow  tube  experiments  and  concludes 
that  laminar  instabilities  will  occur  at  high  heating  levels. 

The  effects  of  free-stream  particulates  on  transition  location  are  not  at 
present  very  well  understood,  particularly  for  the  very  high  Reynolds  numbers 
associated  with  heated  laminar  flow.  Ladd  and  Hendricks  (1982)  have  initiated 
a  quantification  of  these  effects  through  tests  on  a  small  (50  mm  diameter), 
9:1  fineness  ratio  ellipsoid  in  a  water  tunnel.  Their  transition  Reynolds 
numbers  vary  from  3.0  to  7.5  million  and  the  tunnel  water  is  seeded  with  known 
concentrations  of  particles  ranging  from  12  to  133  pm  in  diameter.  Overheats 
range  from  0  to  15°C.  At  the  higher  transition  Reynolds  numbers  considered, 
Ladd  and  Hendricks  conclude  that  the  80  and  133  ym  particles  have  a  degrading 
effect  on  transition  Reynolds  number.  They  estimate  a  critical  particle 
Reynolds  number  (based  on  free-stream  velocity  and  particle  diameter)  to  be  of 
order  700.  It  is  pointed  out  that  this  critical  Reynolds  number  does  depend 
on  heating  level,  surface  finish,  particle  concentration  and  a  judgment  as  to 
where  the  effect  is  important.  Based  on  these  results,  it  appears  that 
free-stream  particles  have  a  detrimental  effect  on  transition  Reynolds  number, 
but  further  work  is  necessary  to  identify  the  mechanisms  by  which  particles 
destabilize  the  laminar  layer.  Chen,  Goland,  and  Reshotko  (1979)  give  a 
nualitative  analvsis  of  particle/boundary  layer  interaction,  but  the  method 
does  involve  some  empiricism.  The  authors  remark  that  detailed  experimental 
efforts  are  needed  to  utilize  the  analvsis  more  fully. 
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The  present  investigation  addresses  the  stability  of  the  laminar  layer 
over  a  relatively  large,  heated  surface  underwater  body.  This  body  has  been 
designed  specifically  for  transition  research  in  the  Garfield  Thomas 
l.22m  diameter  Water  Tunnel  located  at  the  Applied  Research  Laboratory  on 
The  Pennsylvania  State  University  campus.  Considerable  engineering  effort  has 
been  given  to  the  design  and  manufacture  of  this  body  so  that  it  can  be  a 
versatile  tooL  for  transition  research.  The  body  surface  is  heated 
eLectrically  so  that  heating  distributions  can  be  changed  within  rather  broad 
limits.  The  shape  of  the  body  is  one  that  yields  a  mildly  favorable 
pressure  gradient  along  its  length.  The  surface  of  this  body  is  polished 
nickel  and  is  typical  of  "best  effort  '  polishing  by  a  skilled  craftsman. 

The  water  tunnel  in  which  the  heated  body  is  tested  is  adequate  for  many 
laminar  flow  experimental  programs.  Turbulence  management  includes  both 
coarse  and  fine  mesh  honeycomb  and  a  9:1  contraction  ratio  settling  section. 
The  background  noise  is  quite  low  due  to  excellent  cavitation  performance  of 
the  impeller.  The  important  feature  of  this  facility  is  its  high  unit 
Revnolds  number  capability.  After  heating  of  the  water  (48°C),  unit  Reynolds 
numbers  of  34.9  milliou/m  can  be  achieved  at  a  maximum  water  velocity  of 
19.8  m/s. 

The  transition  data  reported  in  this  paper  represent  baseline  data  for  the 
subject  bodv.  One  objective  is  to  estimate,  theoretically,  heat  distributions 
for  the  bodv  based  on  the  analysis  of  Eisenhuth  and  Hoffman  (1981)  that  should 
result  in  a  stable  boundary  laver.  These  estimates  are  thence  compared  with 
the  experimental  findines.  \r.other  objective  is  to  show  qualitatively  any 
dependence  of  transition  Revnolds  number  on  the  concentration  of  free-stream 


particulates  that  are  present  in  the  water 
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2.1  The  1.22m  Diameter  Water  Tunnel 

The  facility  used  in  the  current  heated  body  investigation  is  a  closed 
circuit,  low  turbulence  water  tunnel  which  has  a  test  section  length  of  4.27m 
and  a  diameter  of  1.22m.  A  complete  description  of  this  facility  is  given  by 
Lehman  (1959). 

The  turbulence  in  the  test  section  is  quite  low  for  a  water  tunnel  of  this 
size.  The  settling  section  contraction  ratio  is  9:1  and  its  upstream  section 
contains  a  flow  straightener  and  honevcomb.  The  flow  straightener  is  a  0.61m 
long  square  mesh  arrangement  with  a  cell  size  of  10.16  cm.  The  stainless 
steel  honeycomb  is  situated  0.61m  downstream  from  the  straightener,  is  20.32 
cm  thick,  and  has  an  octagonal  cell  shape  0.56  cm  from  flat  to  flat.  Lumley 
and  McMahon  (1967)  describe  the  design  considerations  for  this  turbulence 
management  system.  The  longitudinal-component  of  turbulence  intensity  in  the 
test  section  and  settling  section  has  been  measured  and  reported  by  Robbins 
(1978).  Figure  1  shows  the  turbulence  intensity  as  a  function  of  test  section 
velocity.  The  higher  levels  below  8  m/s  are  suspected  of  being  due  to  laminar 
instability  within  the  cores  of  the  honeycomb.  The  majority  of  the  data 
presented  in  Section  4  are  for  velocities  greater  than  8  m/s  where  the 
turbulence  intensity  is  of  order  0.1  percent  and  the  average  integral  length 
scale  is  1.25  cm.  As  shown  by  Lauchle,  et  al  (1980),  this  level  is 
sufficiently  low  for  transition  research. 

A  feature  of  the  1.22m  diameter  water  tunnel  worth  noting  is  its  degassing 
and  fiLterir.g  capability.  A  fluid  bypass  svstem  is  available  which  is  used  to 
pass  the  water  through  10  urn  filters  and  de-aeration  modules.  The  maximum 
flow  of  this  svstem  is  190  l/s.  The  de-aeration  modules  are  small,  plastic 


saddle-shaned  units.  Manv  thousands  of  these  modules  are  contained  in  a  large 
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settling  tank  through  which  the  water  flows.  The  resulting  trickle  of  water 
over  the  modules  effectively  removes  dissolved  and  free  gas.  Particulate 
matter  is  removed  by  filters  in  this  bypass  system. 

The  water  temperature  increases  with  time  as  the  facility  is  operated. 
Typically,  for  a  test  section  velocity  of  13  m/s,  the  temperature  increases 
at  a  rate  of  1°C  per  hour.  Testing  at  constant  Reynolds  number  can  be 
achieved  by  computer/tunnel  operator  interaction  because  water  temperature, 
velocity,  and  pressure  data  are  passed  to  an  on-line  VAX  computer  every  10s. 

The  water  quality  can  be  determined  at  any  time  during  the  course  of  a 
test.  The  dissolved  gas  content,  aQ,  is  measured  with  a  Van  Slyke  apparatus. 
This  instrument  is  primarily  used  in  the  medical  diagnosis  of  blood,  but  can 
also  be  used  for  water.  Th^  net  output  is  the  number  of  moles  of  dissolved 
gas  per  total  number  of  moles,  expressed  in  parts  per  million  (ppm). 

The  free  gas  in  the  test  section  can  be  measured  with  a  laser  light 
scattering  system  developed  by  Davis  and  Billet  (1983).  This  system  measures 
the  number  of  bubbles  per  unit  volume  for  each  of  seven  different  resolution 
channels  ranging  from  9  to  48  pm  bubble  diameter. 

Water  purity  is  determined  bv  either  of  two  methods.  The  one  method, 
developed  by  Billet  and  Gates  (1981),  uses  optical  holography  to  sample 
statistically  a  volume  of  the  tunnel  water.  Particle  size  distributions 
covering  the  10  to  70  pm  mean  diameter  range  can  be  determined  with  this 
svstem.  The  other  method  involves  the  passing  of  any  selected  volume  of 
tunnel  water  through  a  N’uclepore®  membrane  filter.  Photomicrographs  of  the 
filter  are  produced  using  a  Scanning  Electron  Microscope  (S.E.M.)  located  at 
the  Materials  Research  Laboratory  on  The  Pennsylvania  State  University  campus. 
The  photomicrograph  is  further  scanned  with  the  S.E.M.  and  a  statistical 
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estimate  of  the  size  distribution  of  particles  is  performed.  This  system  is 
capable  of  counting  any  size  particle,  but  typically  the  measurements  cover  a 
1  to  1000  pm  mean  diameter  range. 

2.2  The  Heated  Body 

The  body  shape  considered  for  this  investigation  is  described 
mathematically  by  a  family  of  modified  ellipse  functions.  It  is  typical  of 
the  kind  of  underwater  body  whose  transition  location  might  be  effectively 
delayed  through  use  of  surface  heating.  The  equation  which  describes  the 
body  coordinates  is  given  by: 

G  x 

y  =  [x(2  -  x)]1^  -  K  — exp  {-  x^/2a^  -  kK  }  +  x^e.  ,  (1) 

n  2az  n 

where  x  =  x/ 1  ,  l  =  nose  length,  y  =  r,  /D  /2,  D  =  maximum  diameter, 
n  n  b  max  max 

K  =  l  /D  ,  a  =  0.3,  Cn  =  0.0303,  and  k  =  0.45227.  Here,  x  is  the  axial 
n  n  max’  ’  u  ’ 

coordinate  and  r^  is  the  body  radius.  The  parameter 

K  C 

=  — n-2°'  exp  {-  l/2a  -  kK.^}  .  (2) 

2a 

Equation  (1)  describes  the  shape  from  x=0tox=£n(£n=  2.44m),  where 
rb  =  0max/2  (Dmax  =  0.32m).  From  x  =  2n  to  x  =  L  (L  =  total  length  =  3.05m) 
the  body  fairs  into  a  0.2m  diameter  sting.  The  fairing  curve  is  a  simple  arc 
of  radius  1.29m  combined  with  an  inflection  curve  at  the  sting  juncture. 
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Calculations  and  measurements  of  the  pressure  distribution  over  this  body 
are  given  by  Lauchle  (1979).  The  computed  potential  flow  results  are  repeated 
here  in  Figure  2.  The  calculations  for  the  pressure  coefficient,  Cp?  ancj 
pressure  gradient  parameter,  8,  account  for  the  solid  blockage  that  occurs  on 
bodies  operating  within  the  confines  of  tunnel  walls.  The  parameter  8  is 
calculated  in  the  standard  way,  i.e.. 


8  = 


2M 

M  +  X  +  1 


(3) 


where 


M  =  (x/u  )(du  /dx) 
e  e 


(4) 


u^  =  edge  velocity  , 

and 

X  =  (2x/r  )(dr  /dx)  .  (5) 

b  b 

It  is  seen  from  Figure  2  that  this  body  supports  a  mildly  favorable  pressure 
gradient  along  most  of  its  length.  The  parameter  8  remains  relatively 
constant  and  positive,  indicating  that  the  pressure  gradient  is  relatively 
constant  and  favorable. 

Given  the  information  contained  in  Figure  2,  temperature  and  heat  flux 
distributions  can  be  predicted  that  should  result  in  laminar  flow  over  the 
body  to  the  maximum  diameter  point  (x  =  1).  These  calculations  are  performed 
using  the  method  of  Eisenhuth  and  Hoffman  (1981).  (Note:  The  numerical 
examples  given  in  this  reference  are  for  this  particular  body.)  The  method 
makes  use  of  similarity  solutions  of  the  boundarv  laver  equations  which 
include  the  heat  transfer  terms,  although  transverse  curvature  terms  are 
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neglected.  The  parametric  procedure  permits  one  to  estimate  heat 
distributions  based  on  either  of  two  criteria.  The  first  criterion  insures 
the  maintenance  of  laminar  flow  by  the  addition  of  just  enough  heat  to  keep 
the  local  displacement  thickness  Reynolds  number  equal  to  the  minimum  critical 
Reynolds  number  established  previously  by  Lowell  and  Reshotko  (1974).  This 
heating  level  is  referred  to  as  the  ''minimum  distribution".  The  second 
criterion  is  based  on  the  fact  that  there  is  a  peak  critical  Reynolds  number 
at  about  65°C  overheat  (Lowell  and  Reshotko  1974).  Thus,  the  method  provides 
enough  heat  so  that  the  peak  critical  Reynolds  number  is  always  maintained. 
This  level  of  heating  is  referred  to  as  the  "maximum  distribution".  Heating 
above  this  level  appears  to  provide  less  enhancement  of  the  laminar  layer 
stability.  The  maximum  temperature  distribution  is  independent  of  velocity 
while  the  minimum  distribution  does  change  with  the  operating  velocity. 

Because  the  local  Reynolds  numbers  are  equal  to  or  less  than  a  critical  value, 
no  amplification  of  spatial  disturbances  in  the  laminar  layer  would  be 
expected.  The  results  of  these  predictions  will  be  discussed  in  Section  4. 

The  heating  system  for  the  test  body  is  designed  to  provide  enough  heat  so 
that  a  theoretical  maximum  heat  flux  distribution  can  be  achieved  with  the 
tunnel  operating  at  its  maximum  velocity  of  19.8  m/s  and  typically  high  water 
temperature  of  32.2°C.  In  order  to  provide  flexibility  in  the  setting  of 
given  heating  distributions,  an  electric  heat  exchanger  is  used  in  the  body. 
This  heat  exchanger  consists  of  thirtv-two  (32)  individually  controlled  axial 
heat  zones.  Figure  3  shows  a  cut-awav  view  of  the  bodv  construction.  A 
typical  heat  zone  is  an  aluminum  ring,  7.b2  cm  in  length  and  2  cm  thick.  The 
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ring  is  fitted  with  a  dense  array  of  commericaliy  available  cartridge  heaters. 
The  array  of  heaters  are  parallel  wired  and  the  common  leads  are  connected 
to  the  power  source  and  controller. 

Each  aluminum  ring,  or  heat  zone,  is  machine  fit  to  its  adjacent  ring. 

When  all  rings  are  stacked  together  (several  internal  tie  rods  secure  the 
assembly)  an  inner  body  of  length  £n  is  formed.  At  first  assembly,  this  inner 
body  was  machine  finished  as  a  unit.  It  was  then  used  as  a  mandrel  for 
electroforming  a  solid  nickel  outer  shell.  The  electrof ormed  nickel  outer 
shell  is  25.4  mm  thick  and  2.44m  long.  Its  surface  is  machined  to  Equation 
(1)  using  a  numerically-controlled  lathe  and  then  buffed  and  polished  to  a 
fine  finish.  Surface  waviness  is  less  than  25  pm  in  10  cm  and  the  rms  value 
of  the  random  surface  roughness  is  of  order  0.5  pm. 

The  inner  rings  can  be  removed  from  the  nickel  shell  for  wiring  or  repair. 
When  the  body  is  reassembled,  a  thermally  conductive,  silver  impregnated 
grease  is  used  between  the  inner  body  and  outer  shell.  Register  marks  permit 
the  body  to  be  reassembled  in  exactly  the  same  way  it  was  during  the 
electroforming  process;  thus,  minimizing  misalignment  which  would  degrade 
heat  transfer. 

The  heated  nose  portion  of  the  body  (0  <  x  <  £n)  is  secured  to  a  stainless 
steel  afterbody  and  sting  mount.  Figure  4  shows  a  photograph  of  the  body  as 
installed  in  the  1.22m  diameter  test  section.  The  afterbody  is  equipped  with 
a  strain  gaged  force  balance  for  total  drag  measurements.  Because  of  certain 
unresolved  problems  with  the  balance,  there  will  be  no  further  discussion  of 


drag  measurements  in  this  paper. 


2.3  Instrumentation 
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The  instrumentation  used  in  the  experimental  program  includes  the 
equipment  necessary  to  measure  body  heat  flux,  body  temperature,  transition, 
water  quality,  and  operational  parameters.  Because  the  techniques  for 
measuring  water  quality  have  been  discussed  in  Section  2.1,  no  further 
discussion  is  given  here. 

The  operational  parameters  include  test  section  velocity,  uQ,  ambient 
pressure,  P0>  and  ambient  water  temperature,  T0.  Velocity  is  determined  from 
the  pressure  drop  between  the  settling  section  and  the  entrance  of  the  test 
section.  The  ambient  pressure  is  measured  at  the  test  section  entrance. 

These  pressures  are  determined  using  Validyne  AP-10  absolute  and  DP-15 
differential  pressure  transducers.  The  water  temperature  is  monitored  in  the 
settling  section  with  an  RTD  sensor  with  digital  output.  All  signals  are 
passed  through  a  100  channel  Data  Systems  multiplexer  and  then  to  a  VAX 
digital  computer  system  for  subsequent  conversion  to  engineering  units.  The 
multiplexer  updates  all  readings  every  10s. 

As  noted  in  Figure  3,  each  heat  zone  contains  a  group  of  thermocouples 
(typically  four).  These  thermocouples  are  made  from  chromel/alumel ,  ANSI  Type 
K  wires.  The  junction  is  crimped  using  commercial  nickel  sleeves.  The 
measurement  junction  is  very  close  to  the  outside  surface  of  the  inner  shell. 

A  procedure  for  estimating  wall  temperature  from  measurements  made  at  the 
inner/outer  shelL  interface  is  discussed  in  Section  3.  The  ninety  (90)  body 
thermocouo Les  pass  from  the  bodv  to  a  copper  zone  box  along  with  a 
thermocouple  from  an  Omeea  Model  TRC-III  ice  point  cell  and  then  to  the  100 
channel  multiplexer.  The  voltages  are  monitored  with  a  H/P  voltmeter  and 
transmitted  to  the  VAX  computer.  The  wires  are  individually  calibrated  prior 
to  installation,  where  the  conversion  of  voltage  to  temperature  units  is 
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performed  In  the  VAX  computer.  CRT  display  of  the  body  temperature 
distribution  is  achieved  in  nearly  real  time.  The  accuracy  of  the  temperature 
measurement  system  is  estimated  to  be  ±  0.5°C. 

Body  heat  flux  distributions  are  estimated  from  the  ratio  of  the 
electrical  power  supplied  to  each  of  the  thirty-two  (32)  heat  zones  to  the 
wetted  surface  area  of  the  given  zone.  This  procedure  is  most  valid  under  the 
conditions  of  forced  convection  when  the  heat  flux  vector  is  radially  outward. 
It  has  been  estimated  (Stinebring  1977)  that  for  the  typical  water  velocities 
used,  the  ratio  of  outward  heat  flow  to  inward  (toward  center  of  bodv')  heat 
flow  is  of  order  10  .  The  electrical  power  of  each  zone  is  varied  with 
custom-designed  SCR-type  controllers;  each  rated  at  7.3  kW.  The  power  can  be 
adjusted  in  one  (1)  percent  increments  of  the  load. 

A  Magtrol  digital  power  meter  has  been  modified  by  General  Electric  Corp. 
to  selectively  monitor  the  power  supplied  to  each  heat  zone.  This  device 
integrates  the  product  of  voltage  and  amperage  for  a  period  of  30s  and 
displays  the  wattage  digitally.  The  manufacturer  of  the  power  meter  states 
the  accuracy  to  be  ±  0.25  percent  of  the  displayed  reading.  For  the  typical 
power  levels  used  in  this  investigation,  an  estimated  accuracy  is  of  order 
±  6  W. 

In  order  to  maintain  a  high-quality  surface  finish,  the  hot-film  probes 
used  for  monitoring  transition  are  not  placed  in  the  heated  portions  of  the 
bodv.  Rather,  three  flush  mounted  films  (TSI  Model  1471)  are  spaced  at 
120°  increments  around  the  body  at  the  downstream  end  of  the  last  heat  zone. 
The  arc  length  distance,  s,  to  the  probes  is  2.12m.  Conseauent ly ,  the 
transition  Reynolds  number  is  measured  for  a  fixed  distance  and  is  varied  bv 
varving  the  free-stream  velocity.  The  hot-film  signals  are  processed  using 
standard  DISA  55MOI  anemometers  and  DTSA  55026  signal  conditioners. 
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Transitional  flow  over  a  fixed  hot-film  probe  is  intermittent  due  to  the 
creation  and  convection  of  turbulent  spots.  This  intermittent  flow  is 
described,  statistically,  by  the  intermittency  factor,  y,  and  the  burst 
formation  rate,  N.  An  Indicator  function,  I(t),  is  defined  as  a  zero-one 
function;  it  is  zero  if  the  boundary  layer  is  laminar  and  it  is  one  if  the 
layer  is  turbulent.  The  factor  y  Is  the  time  average  va^ue  of  I(t),  and  N  is 
a  count  of  the  average  number  of  step  changes  in  I ( t )  per  unit  time.  An 
intermittency  detector  is  used  to  generate  I(t).  Its  operation  is  as  follows 
(l)  The  hot-film  signal  is  high-pass  filtered  to  eliminate  very  low  frequency 
components  (typically  set  at  50  Hz),  (2)  the  signal  is  differentiated  and 
full-wave  rectified  so  as  to  enhance  high-frequency  turbulent  fluctuations, 
and  be  representative  of  a  fluctuating  vorticity  signal  (Taylor's  frozen 
vorticity  assumption),  (3)  it  is  integrated  with  a  short  time  constant 
low-pass  filter  in  order  to  smooth  the  fluctuations  within  a  burst,  and  (4) 
the  conditioned  signal  then  passes  through  a  Schmidt  trigger  with  adjustable 
threshold.  The  output  of  the  Schmidt  trigger  is  I(t).  Figure  5  shows  these 
steps  as  performed  on  a  typical  hot-film  signal  obtained  in  a  transitioning 
boundary  layer.  The  intermittency  factor  is  determined  by  passing  I(t) 
through  an  integrating  voltmeter  and  N  is  established  by  passing  I(t)  through 
a  counter  that  triggers  only  on  the  leading  edge  of  the  unit  step  functions. 
Both  quantities  are  determined  for  record  lengths  of  100s  or  more.  As 
explained  below,  the  value  of  y  is  used  to  establish  transition  for  these 
experiments . 

3.  EXPERIMENTAL  PROOF, PORES 

The  hot-film  probes  are  used  onlv  for  the  detection  of  transition. 
Therefore,  an  absolute  calibration  is  not  necessarv.  The  intermittency 
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detector  is  set  up  using  a  storage  oscilloscope.  The  velocity  is  adjusted 
so  that  an  approximately  50  percent  intermittent  signal  is  observed.  The 
detector  threshold  voltage  is  adjusted  until  the  indicator  function  tracks 
with  the  turbulence  bursts  sensed  by  the  films.  This  setting  is  very  rareiv 
changed  through  the  course  of  a  test.  This  stability  is  attributed  to  the 
fact  that  the  DC-component  of  the  film  signal  is  removed  by  high-pass 
filtering. 

Transition  Reynolds  number  is  determined  from  measurements  of  y  vs.  uQ. 

The  velocity  is  set  at  a  low  value  where  y  =  0  (fully  laminar  flow  over  the 
hot  film  probe).  It  is  then  increased  incrementally  and  fixed  while  y  and  N 
are  established.  After  the  higher  velocities  (where  y  =  1)  have  been 
achieved,  the  sequence  is  repeated  by  decreasing  uQ  in  increments.  The  values 
of  y  and  N  are  also  determined  during  the  decreasing  velocity  part  of  the  run. 
It  is  again  noted  that  100s  (or  more)  of  integration  time  is  used  for  each 
datum.  The  abscissa  is  converted  to  arc  length  Reynolds  number  (Res  =  uQs/v, 
where  v  is  the  kinematic  viscosity  determined  for  the  water  temperature  at  the 
time  the  datum  is  measured,  and  s  =  2.12m).  The  arc  length  Reynolds  number  at 
the  point  where  y  =  50  percent  is  defined  as  the  transition  Reynolds  number. 
These  numbers  may  be  multiplied  by  1.085  to  obtain  the  arc  length  Reynolds 
number  based  on  edge  velocity,  ug.  This  constant  arises  from  the  value  of 
Cp  at  s  =  2.12m. 

The  heat  f Lux  distributions  are  measured  in  a  straightforward  manner.  As 
previously  noted,  the  electrical  power  supplied  to  a  given  heat  zone  is  simply 
divided  bv  the  wetted  surface  area  of  the  zone.  The  procedure  is  approximate 
because  of  the  small,  unknown  amount  of  heat  that  leaks  axial Iv  from  zone  to 
zone.  N'o  attempt  has  been  made  to  account  for  this  effect. 
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The  surface,  or  wall  temperature,  Tw,  is  estimated  from  the  measured 
temperature,  Tm,  by  an  elementary  method.  Figure  6  shows  a  cross  section  of 
the  body  in  a  plane  containing  a  thermocouple.  If  q  is  the  heat  flux  for  the 
particular  zone,  and  £  is  the  axial  Length  of  the  zone, 

T  =  T  -  —3—  [£n(r  /r  )/k  +  £n(r  /r  )/k  ]  .  (6) 

w  m  2it£  2m  AL  w  2  NI 

The  constants  and  kfjj  are  the  thermal  conductivities  of  aluminum  and 
nickel,  respectively.  Equation  (6)  assumes  a  perfect  contact  at  the 
aluminum/nickel  interface  and  that  the  heat  flux  is  predominantly  radially 
outward. 

Free  bubble  distributions  are  measured  with  the  light  scattering  nucieii 
counter  of  Davis  and  Billet  (1983).  These  measurements  are  performed  in  the 
test  section  near  the  maximum  diameter  point  of  the  body.  Because  bubble 
diameter  is  affected  by  dynamic  pressure,  these  distributions  are  measured 
over  a  range  of  velocities.  The  tunnel  pressure,  PQ,  has  been  held  constant 
at  241.3  kPa  (absolute)  for  all  data  presented  in  this  paper. 

The  water  samples  used  in  determining  particulate  distributions  are 
collected  through  a  valve  located  on  the  bottom  side  of  the  sett  Ling  section. 
Lauchle  and  Crust  (1980)  report  that  the  measured  distributions  depend  very 
little  on  free-stream  velocity  or  pressure. 

4.  RESULTS 

4 . 1  Heat  Distributions 

The  heat  distributions  selected  for  experimental  studv  are  those  minimum 
distributions  calculated  for  uQ  =  3.05,  7. 82,  12.19,  and  16.78  m/s  at  a  water 
temperature,  TQ,  of  25.6°C.  As  noted  earlier,  each  of  these  distributions 
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should  provide  enough  heat  to  keep  the  boundary  layer  neutrally  stable  to  the 
ooint  of  maximum  diameter  on  the  body  at  velocities  less  than  .or  equal  to  the 
value  quoted.  The  heat  flux  distributions  are  shown  in  Figure  7  while  the 
corresponding  temperature  distributions  are  given  in  Figure  8.  Integration  of 
the  heat  flux  distribution  over  the  wetted  surface  area  of  the  body  gives  the 
total  heating  power  supplied  to  the  body.  These  integrated  values  are  denoted 
by  Qj  and  are  indicated  on  Figures  7  and  8.  In  later  discussions  of  the 
transition  data,  Q-p  will  be  used  as  the  parameter  that  identifies  which 
distribution  is  being  referred  to.  It  is  seen  from  Figure  7  that  the  measured 
heat  flux  distributions  agree  well  with  the  theoretical  distributions. 

The  agreement  between  the  measured  wall  temperature  and  that  predicted 
is  fair  (Figure  8).  The  temperature  data  shown  are  measured  at  the  value  of 
u 0  noted  and  under  laminar  flow  conditions.  These  data  include  the  correction 
for  skin  thickness  described  by  Equation  (6),  in  which  measured  values  of  q 
are  used.  The  scatter  in  these  data  is  attributable  to  the  slight  scatter  in 
the  heat  flux  distributions  and  to  the  non-perfect  contact  between  the  inner 
and  outer  shells  of  the  body. 

Temperature  data  are  not  shown  for  the  two  highest  heating  conditions 
because  the  boundary  layer  was  observed  to  be  intermittent  at  the  hot-film 
location.  This  intermittency  greatly  affects  the  heat  transfer 
characteristics  of  the  body  and  results  in  random  fluctuations  of  the  wall 
temperature.  Reoeatable  data  couLd  not  be  produced  near  the  end  of  the  body 
where  the  Laver  is  intermittent.  Although  not  explored  in  detail,  the 
thermocouple  data  can  be  used  to  indicate  transition  Location.  However,  the 
axial  resolution  is  poor  when  the  intermittent  flow  regimes  that  precede 
f u L Lv-deve Loped  turbulent  flow  are  Large. 
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4.2  Water  Quality 

One  of  the  objectives  of  this  studv  is  to  establish  a  potential 
correlation  between  the  amount  of  entrained  particulate  matter  in  the  water 
and  the  transition  Reynolds  number  of  the  body  under  various  levels  of  surface 
heating.  An  effective  wav  of  doing  this  is  to  start  out  with  relatively  dirty 
water  and  then  obtain  transition  data  at  selected  intervals  during  a  filtering 
process.  The  tunnel  water  (400  m^  total  volume)  is  drained  and  replaced  with 
fresh  water  from  a  local  well.  This  water  is  traditionally  high  in 
particulate  matter  and  dissolved  gas.  Through  use  of  the  tunnel  bypass 
system,  the  particles  and  dissolved  gas  are  removed,  but  the  process  requires 
several  hours  to  significantly  reduce  the  levels.  Furthermore,  the  bypass 
operation  can  be  stopped  at  any  time  and  the  facility  operated  in  its  routine 
manner  while  transition  data  are  collected. 

A  history  of  the  water  quality,  as  measured  in  this  program,  is  given  in 
Figure  9.  Here,  total  concentration  is  defined  as  the  total  number  of 
particles  or  bubbles  measured  without  regard  to  size,  per  unit  volume.  Open 
circLes  denote  those  particles  measured  using  the  holographic  technique  and 
the  solid  circles  represent  additional  data  acquired  using  the  S.E.M. 
technique.  Triangles  are  used  to  denote  the  free  gas  bubbles  measured  with 
the  Light  scattering  method.  The  mean  diameter  of  the  particLes  and  bubbles 
is  also  shown  in  Figure  9  as  a  function  of  the  filtering  time.  Transition 
data  are  obtained  for  each  condition  indicated  by  a  datum  point  on  this 
figure.  Tt  is  emphasized  that  a  given  condition  was  maintained  for  as  long  a 
period  as  necessary  in  order  to  acquire  the  transition  data. 

The  purest  water  considered  occurs  after  several  hours  of  filtering. 

The  concentration  of  bubbles  and  particles  are  nearly  equal  at  this  condition, 
biit  very  Low;  less  than  lO/cm^.  Prior  to  this  condition.  Figure  9  shows 
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that  the  water  Is  dominated  by  particLes,  typically  20  to  26  ym  in 
diameter. 

Measured  particLe  and  bubble  distributions  are  shown  in  Figure  10.  The 
cross-banded  area  of  this  figure  represents  the  spread  of  typicaL  particulate 
distributions  measured  in  various  oceans  (Barker  and  Gile  1981).  From  this 
figure  it  is  seen  that  the  water  tunnel  distributions  considered  fall  below 
and  within  those  of  typical  ocean  water.  Error  bands  are  indicated  tor  the 
bubble  distributions  because  these  measurements  exhibit  considerable  scatter. 
The  scatter  is  due  to  the  dependence  of  bubble  size  on  velocity. 

4 . 3  BaseLine  Transition  Results 

In  this  section,  transition  data  are  presented  for  the  purest  conditions 
of  tunnel  water  (t/tT0TAL  >  0.94)  where  the  aggregate  concentration  of 
particLes  and  bubbles  is  Less  than  15/cm^.  These  data  are  considered  baseline 
data  and  represent  the  best  performance  of  the  body.  It  is  noted  that  the 
body  was  installed  in  the  tunnel  four  different  times  during  a  12  month  time 
period,  and  that  these  baseline  transition  data  were  repeatable  to  within 
5  percent.  The  data  presented  are  for  one  hot-film  probe  only  because  two  of 
the  three  probes  installed  in  the  body  did  not  perform  satisfactorily. 

Figures  l l  through  14  show  the  baseline  transition  region  statistical 
quantities  y  and  N  as  functions  of  the  arc  length  Reynolds  number.  The 
transition  Reynolds  number  is  identified  at  7  =  50  percent.  The  intermit cency 
and  burst  rate  distributions  are  seen  to  broaden  as  the  heating  level  (and 
transition  RevnoLds  number)  increases.  Within  the  experimental  scatter, 
no  dependence  of  y  (or  N)  on  whether  u0  is  increasing  or  decreasing  is 
observed.  A  hysteresis  might  be  expected,  however,  if  the  integration  time 
were  decreased  and  the  rate  of  velocity  change  increased.  Such  data  could  no 
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longer  be  classed  as  steady  state,  but  as  transient  and  would  depend  on  the 
thermal  inertia  of  the  body  itself. 

The  cold  body  intermittency  distribution  is  considered  to  be  quasi_steady . 
With  the  hot-film  orobe  situated  at  a  relatively  large  arc  length  distance 
from  the  nose,  the  velocity  must  be  set  quite  low  (uQ  ~  1.5  m/s)  for  laminar 
flow.  At  speeds  Lower  than  3  m/s,  the  water  velocity  drifts  thus 
necessitating  the  collection  of  data  over  shorter  time  intervals. 

Nevertheless,  the  cold  body  transition  Reynolds  number  of  4.5  million  agrees 
weil  with  that  measured  by  other  methods  on  another  body  of  this  size  and 
shape  (Lauchle,  et  al  1980). 

The  heating  distributions  associated  with  the  transition  data  of 
Figures  11  through  14  are  based  on  predictions  that  require  the  local 
displacement  thickness  Reynolds  number  to  be  less  than  a  critical  value, 
R'e6*crit’  A  8i-ven  heat  distribution  should  therefore  provide  enough  heat 
to  maintain  full  laminar  flow  over  the  body  at  the  velocity  noted  for  the 
prediction.  A  critical  arc  length  Reynolds  number  can  be  defined  as  one  based 
on  the  velocity  used  in  the  heat  prediction  and  the  arc  length  from  the  nose 
to  the  hot-fiLm  probe.  A  comparison  of  this  critical  Reynolds  number  with  the 
measured  transition  Reynolds  number  is  given  in  Figure  15.  One  would  expect 
the  critical  curve  to  fall  substantially  below  the  transition  curve,  as  it 
does  for  the  lower  values  of  Qq-.  However,  for  Ox  >  60  kW,  the  two  curves 
converge  and  cross.  This  behavior  is  not  inconsistent  with  what  Barker  and 
Gi  le  (1981)  observed  for  a  heated  pipe  flow.  The  theoretical  methods  used  by 
Barker  and  Gi  le  are  analogous  to  those  used  here  (F.isenhuth  and  Hoffman  1981). 
It  is  reasoned  that  for  the  higher  Reynolds  numbers  (or  higher  heating 
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levels),  free-stream  disturbances,  surface  flaws,  and  system  idiosyncrasies 
become  more  important;  their  presence  affects  the  measured  results,  but  are 
unaccounted  for  in  the  theoretical  models. 

The  application  of  Emmons'  (1951)  spot  theory  allows  one  to  express  the 
burst  frequency  in  terms  of  the  intermittency  factor.  Following  Gedney 
(1979),  a  line  source  density  function  is  assumed,  for  which 

NAx/uq  =  4.76  [(1  -  YUn(^)]1/2  .  (7) 

The  streamwise  extent  of  the  transition  zone  is  denoted  by  Ax.  This  length 
is  not  known  for  the  present  measurements,  so  it  must  be  estimated.  Chen 
and  Thyson  (1971)  suggest  that 

2/3 

ReA  =  60Re .  ,  (8) 

Ax  t  ’ 


where  Ret  is  the  length  Reynolds  number  where  bursts  first  appear. 
Rearrangement  of  Equation  (7)  using  Equation  (8),  and  using  the  fact  that 
Res  ~  Ret  (s  =  distance  to  the  measurement  probe)  permits  one  to  write: 


Re 


Ns 

17T 


-  0t[  (  1  - 


1/2 


y)£n(yr^)] 


(9) 


s  o 

where  a  is  a  proportionality  constant. 

Figure  16  shows  the  baseline  non-dimensional  burst  rates  as  a  function  of 
intermittency  factor  for  the  four  heating  levels  considered.  The  constant 
a  of  Equation  (9)  is  adjusted  so  that  the  theoretical  curve  matches  the 
data.  This  matching  reveals  that  at  depends  on  Qj  in  a  non-linear  manner. 
Figure  17.  This  dependence  cannot  be  exDlained  further,  until  additional 
information  is  obtained  on  the  behavior  of  Ax  with  varving  levels  of 
svirface  heating.  (The  constants  ir.  F.qs.  (7)  and  (8)  are  based  on  experiments 


with  unheated  waiis.) 
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4 .4  Effect  of  Increased  Particulate  Level 

In  reference  to  Figure  9,  the  baseline  transition  results  of  Section  4.3 
are  for  normalized  filtering  times  greater  than  0.94.  In  this  section, 
transition  data  are  presented  for  much  shorter  filtering  times  where  the  water 
is  contaminated  with  high  levels  of  particulate  matter.  Because  the  number  of 
particles  per  unit  volume  is  far  greater  than  the  number  of  observed  bubbles, 
the  particle  concentration  can  be  used  as  an  independent  variable. 

Figure  18  shows  the  measured  transition  Reynolds  numbers  for  the  body 
operating  at  three  different  heating  levels  as  a  function  of  total  particle 
concentration.  Here,  total  concentration  is  the  total  number  of  particles 
measured  throughout  the  10  to  70  pm  diameter  size  range  per  unit  volume  of 
water.  Figure  18  shows  that  the  laminar  flow  performance  of  the  body  degrades 
with  increasing  concentration  of  free-stream  particulates.  The  effect  seems 
to  be  progressively  more  important  as  the  heating  level  increases.  This  is 
more  easily  seen  by  replotting  the  data  as  shown  in  Figure  19.  The  transition 
Reynolds  number  is  decreased  by  as  trajch  as  30  percent  due  to  the  introduction 
of  the  higher  levels  of  free-stream  particles. 

The  results  given  in  Figures  18  and  19  show  that  the  decrease  in 
transition  Reynolds  number  with  increasing  particle  concentration  is  a 
gradual  dependence.  This  is  analogous  to  what  Klebanoff,  Schubauer,  and 
Tidstrom  (1955)  observed  for  two-dimensional  roughnesses  on  a  flat  plate. 

The  transition  Reynolds  number  gradually  decreases  as  the  roughnesses 
protrude  farther  and  farther  into  the  boundary  layer.  The  mechanism  is 
explained  simply.  The  roughness  creates  an  inflected  velocity  profile 
which  leads  to  wave  motion  in  the  laminar  layer.  For  small  roughnesses, 
these  waves  damp  out,  but  as  the  roughness  height  increases,  the  waves 
are  less  damped  and  grow.  The  amplification  of  these  roughness-induced 
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waves  leads  to  earlier  transition.  It  is  noted  that  Klebanoff,  et  al  (1955) 
observed  a  different  effect  for  three-dimensional  spherical  elements  that 
were  fixed  on  the  surface.  A  very  pronounced  change  in  transition  Reynolds 
number  occurs  when  the  sphere  diameter  becomes  approximately  equal  to  6*,  the 
displacement  thickness.  From  these  three-dimensional  results,  an  average 
critical  sphere  diameter  Reynolds  number  (based  on  local  velocity  at  the  top 
of  the  sphere)  of  577  is  suggested.  The  mechanism  for  the  abrupt  change  in 
transition  Reynolds  number  appears  in  the  formation  of  an  attached  turbulent 
wedge  from  the  roughness  element.  Thus,  transition  of  the  boundary  layer  is 
caused  by  a  direct  seeding  of  turbulence  from  the  wake  of  the  attached 
element,  rather  than  from  the  amplification  of  wavelike  disturbances  as  were 
observed  to  occur  for  the  two-dimensional  elements.  Hall  (1967)  reports 
similar  results  for  spheres  that  are  attached  as  well  as  detached  from  the 
surface. 

In  light  of  the  above  discussion,  the  current  results  suggest  that  free- 
stream  particles  that  enter  the  boundary  layer  affect  transition  by  inducing  a 
locally  inflected  velocity  profile.  Exactly  where  along  the  body  this  effect 
becomes  important  cannot  be  established  from  the  data  obtained,  but  it  is 
speculated  that  the  particle  size  relative  to  the  displacement  thickness  is  a 
key  parameter.  Calculations  of  the  heated,  laminar  boundary  layer  development 
along  this  body  have  been  performed  using  the  methods  of  Gentry  and  Wazzan 
(1976).  For  the  75.4  kW  heating  condition,  and  for  a  velocity  of  14.6  m/s, 
which  is  approximate Iv  the  observed  baseline  transition  velocity,  6*  is 
calculated  to  be  106  pm  at  s  =  1.92m  (iusl  upstream  of  the  transition  probe). 
At  this  station  along  the  bodv,  the  calculation  suggests  that  dmax/ 6*  -  0.66, 
which  may  not  be  an  insignificant  value.  At  stations  closer  to  the  nose. 
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this  ratio  would  be  expected  to  be  larger,  but  again,  it  cannot  be  determined 
where  along  the  body  a  particle  may  actually  enter  the  boundary  layer.  This 
is  a  point  which  possibly  deserves  further  investigation. 

As  a  final  note,  the  observance  of  no  abrupt  change  in  transition 
Reynolds  number  with  changes  in  particle  concentration  (10  pm  <  d  <  70  pm) 
suggests  that  turbulent  seeding  from  the  particles  does  not  occur  in  the 
experiments  described  here.  Such  seeding  may  possibly  occur,  however,  if  some 
of  the  particles  actually  stick  to  the  surface.  As  pointed  out  by  Talbot, 
et  al  (1980),  thermophoresis  considerations  would  make  this  possibility 
unlikely  for  small  particles  in  a  boundary  layer  over  a  heated  surface.  With 
no  clearly  defined  "critical"  concentration,  a  critical  particle  Reynolds 
number  cannot  be  determined  from  the  results  of  these  experiments. 

5.  CONCLUSIONS 

The  present  experimental  investigation  has  shown  that  the  transition 
Reynolds  number  can  be  significantly  increased  on  a  heated,  underwater  body. 
The  maximum  transition  Reynolds  number  achieved  is  36.4  million  with  an 
average  overheat  of  25°C.  The  total  heating  power  for  this  level  of 
performance  is  93.3  kW.  The  maximum  Reynolds  number  can  be  compared  to  4.3 
million  for  the  body  operating  with  no  heat.  Although  higher  heating  levels 
were  not  considered,  the  results  indicate  that  further  increases  in  heating 
power  would  yield  only  marginal  increases  in  the  transition  Reynolds  number. 
For  heating  levels  below  60  kW  (average  overheat  20°C),  the  increase  in 
transition  Reynolds  number  with  heat  agrees  well  with  theory.  Above  this 
level,  the  agreement  is  not  good.  A  similar  trend  was  observed  for  a  heated 
pipe  flow  (Barker  and  Gile  1981).  This  degradation  in  performance  is 
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attributed  to  free-stream  disturbances,  surface  disturbances,  and  system 
idiosyncrasies,  none  of  which  are  accounted  for  in  the  theory,  and  none  of 
which  could  be  quantitatively  determined  in  the  context  of  these  experiments. 

Particles  in  the  free  stream  are  considered  a  free-stream  disturbance  that 
may  lead  to  increased  laminar  instability.  Their  effect  on  transition 
Reynolds  number  has  been  demonstrated  in  these  experiments.  The  particles 
considered  range  from  10  pm  to  70  pm  in  diameter.  At  high  concentrations 
(~  200/cm^),  these  particles  decrease  the  transition  Reynolds  number  by 
as  much  as  30  percent  at  high  heat  levels.  The  degradation  is  less  severe 
at  lower  heating  levels.  Because  transition  Reynolds  number  increases 
monotonically  with  heating  level,  an  alternative  interpretation  is  that  the 
instabilities  created  by  a  given  size  range  of  particles  becomes  progressively 
more  severe  as  the  particle-free  (baseline)  transition  Reynolds  number 
increases.  This  is  an  expected  trend  because  the  ratio  of  particle  size  to 
laminar  boundary  layer  thickness  becomes  larger  at  a  fixed  arc  length  as  the 
Reynolds  number  increases. 

The  decrease  in  transition  Reynolds  number  with  increasing  concentration 
of  particles  is  a  gradually  decreasing  function.  It  is  difficult,  therefore, 
to  identify  a  critical  concentration  from  which  critical  particle  Reynolds 
numbers  can  be  inferred.  Additional  work  is  necessary  to  identify  these 
particle  Reynolds  numbers.  The  experiments  are  best  performed  at  constant 
Reynolds  number  while  various  sizes  of  particles  are  introduced  at  a  fixed 
concentration.  Work  is  continuing  in  this  area. 
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Figure  4.  Photograph  of  the  heated,  laminar  flow  body  in  the  1.22m  diameter 
test  section. 


Figure  5.  Oscilloscope  traces  showing;  (a)  Hot  film  sig 
after  differentiation  and  full-wave  rectificat 
after  low-pass  filtering,  and  (d)  Schmidt  trig 
(indicator  function)  after  o,perati»n  on  Signal 


measured  heat  flux  distributions  for  the  heated 
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Figure  8.  Predicted  and  measured  wall  temperature  distribui  _>ns  for  the 
heated  body. 
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Variation  of  free-stream  particulates  and  bubbles  with  bypass 
filtering  time. 
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BUBBLE  OR  PARTICLE  DIAMETER,  d  (pm) 


Figure  10.  Particle  and  bubble  size  distributions  measured  during  various 
phases  of  the  experimental  program. 
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Figure  14.  Variation  of  turbulence  intermittency  and  burst  rate  with  Reynolds  number  for  baseline 
conditions  and  Q„  =  93.3  kW. 
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Figure  15 
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Comparison  of  experimentally  determined  transition  Reynolds 
numbers  with  theoretical  critical  Reynolds  numbers. 
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